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Fluctuations in fluid velocity and fluctuations in electric fields may both give rise to forces acting 
on small particles in colloidal suspensions. Such forces in part determine the thermodynamic stability 
of the colloid. At the classical statistical thermodynamic level, the fluid velocity and electric field 
contributions to the forces are comparable in magnitude. When quantum fluctuation effects are 
taken into account, the electric fluctuation induced van der Waals forces dominate those induced 
by purely fluid mechanical motions. The physical principles are applied in detail for the case of 
colloidal particle attraction to the walls of the suspension container and more briefly for the case of 
forces between colloidal particles. 

PACS numbers: 82.70.Dd, 82.30.Dd, 83.80.Hj 



I. INTRODUCTION 

Fluctuations in thermodynamic field parameters may 
give rise to forces acting on colloidal particles. Examples 
of such field parameters include velocity fields within the 
fluid surrounding the particle and local electric fields. 
Our purpose is to examine the magnitudes of these two 
force effects. We seek to ascertain which of these forces 
have the dominant effect on the thermodynamic proper- 
ties of colloidal suspensions. 

The classical statistical thermodynamic contribution 
to fluctuations forces are scaled by the thermal energy 
kgT. The forces due to quantum mechanical zero point 
fluctuations involve a frequency scale Woo of motion and 
are thereby scaled by the energy quantum TkJoo. The 
total magnitude of such forces depend upon whether or 
not the fluctuations are classical or quantum mechanical 
in nature. We find for temperatures in the neighborhood 
of room temperature that the electric field fluctuation 
forcepj dominates the hydrodynamic fluctuation forceQ 
in the fully quantum mechanical theory. 

The Einstein theory of purely classical statistical ther- 
modynamic fluctuation forces will be reviewed in SeclTTl 
The general results are employed for the specific exam- 
ple of the force on a colloidal spherical particle due to 
a neighboring hard wall. The hydrodynamic case is dis- 
cussed in Sec IIIII while the electric field case is discussed 
in Sec llVI The general theory of quantum mechani- 
cal fluctuations are explored in Secfvland the frequency 
scales of both fluid mechanical and electrical fluctuations 
are considered in Sec I VII The frequency scales are such 
that the fluid mechanical fluctuation forces are classi- 
cal while the electric dipole fluctuation forces are quan- 
tum mechanical. The latter thereby dominate the for- 
mer as discussed in Sec lVIII Although the electric field 
static quantum fluctuation forces dominate the classical 
fluid static fluctuation forces, the fluid forces arc nonethe- 
less observable. If the bandwidth of experimental obser- 
vations of Brownian motion coincides with a frequency 
regime wherein fluid mechanics holds true, then fluid fluc- 
tuation forces may be (and have been) measured. This 



point is discussed in Sec lVIIll wherein the formula for 
the fluid and electrical fluctuation forces are exhibited 
for two well separated spheres. In the concluding Sec lIXI 
the electric fluctuation forces are shown to more strongly 
determine the phase properties of a colloid, i.e. whether 
the colloidal particles will form a smooth colloidal sus- 
pension or whether the colloidal particles will undergo 
phase separation. 



II. STATISTICAL THERMODYNAMICS 

Let us consider two sets of coordinates, Q = 
{Qi, . . • , Qm) and X = (Xi, . . . , X n ). For the moment, 
let us fix Q and assume that X undergoes classical ther- 
mal fluctuations 3j. The Einstein probability^ P of ex- 
hibiting a deviation AX = (AAi, . . . , AX n ) from ther- 
mal equilibrium has a Gaussian form determined by an 
activation free energy AF. The activation free energy is 
the minimum isothermal work done on the system by the 
environment in order to to produce the fluctuation AX; 
Whence it follows that 

P cx e - Ai? / feBT , 
1 " 

AF = 2T, G ^k(Q)^AX k , (1) 

wherein the matrix ||Gjfc(Q)|| describes an effective 
Hook's law compliance. The static classical fluctuation- 
response theorem dictates that 



AX i AX j =k B TG ij (Q), (2) 

which follows directly from Eq.QJ. The forces / = 
■ ■ ■ , fm) conjugate to Q = (Qi, . . . , Q m ) may be de- 
rived from the free energy 

3,k 
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Hard Wal 




FIG. 1: Shown is a metal sphere of radius a submerged in a 
fluid at a distance z from a hard wall. We seek to compute the 
fluctuation forces on the sphere due to the presence of the wall 
in part to establish the physically dominant contributions. 

The mean value of this fluctuation force 

may be evaluated via Eq.JSJ as 

fi - -—)^G k] (Q) ^ , 

fl = ^^ln(det||G^(Q)||). (5) 

Employing a reference matrix \\G j k \\ which is indepen- 
dent of Q, we may write Eq.(|SJ) in terms of an effective 
potential U(Q), i.e. 

- = dU(Q) 
h dQi ' 

U(Q) = M^lndet(G- 1 (Q)GW) . (6) 
If the G(Q) obeys the Green's function matrix equation 
G(Q) = G^ +G(°)E(Q)G(Q), (7) 

then 

U(Q) = ^lndet(l-G (0) £(Q)) 

U{Q) = -^f^Tr(G(°)E(Q)) n . (8) 

The central result of this section resides in Eqs.©, @ 
and © which describe the effective potential U(Q) of 
classical statistical thermodynamic fluctuation forces in 
terms of the determinants of the compliance matrices. 



The quantum mechanical version of fluctuation forces 
will be briefly discussed in the following work where we 
will explicitly compute the situation shown in Fig^ We 
consider two contributions to the fluctuation forces be- 
tween the sphere and the wall, namely (i) fluid velocity 
fluctuations and (ii) electric field fluctuations. 



III. FLUID VELOCITY FLUCTUATIONS 

For a single metal sphere (in a fluid) with momentum 
p and position r, the minimum isothermal work required 
to produce the momentum is given by the total kinetic 
energy 

V ' j,k 

where the mass matrix ||Mjfe(r)|| plays the role of the 
compliance in Eq.QJ. The static fluctuation response 
Eq.(|5J now reads as the equipartition theorem 

fjpt = k B TM jk (r). (10) 

Let us consider the mass matrix | \Mjt- (r) 1 1 in more detail. 

If the sphere were very far from the wall z = oo, then 
the mass matrix would be given by 

M ( % = 5 jk (M + M ) = 5 jk (M + ^) , (11) 

wherein p is the mass density of the fluid and M is the 
mass of the sphere. The Euler mass which is half the 
mass of the displaced fluid, has the following physical 
interpretation^ . The kinetic energy of a sphere moving 
slowly with velocity v through an infinite bulk fluid is 
given by 

1 3 

Ktotal = Kpartide + Kfi uid = - ^ M j k VjV k 

j,k 

K total = lAflvp + ^lvl 2 . (12) 

The mass M enters into the particle kinetic energy 
Kpartide = (l/2)^l v | 2 - If the particle moves through 
the fluid, then the fluid exhibits a dipolar back flow con- 
tribution Kfi U id = (1/2)m|v| 2 to the total kinetic energy. 

When the sphere is at distance z < oo from the wall, 
the back flow fluid mass current vector must have a zero 
component normal to the boundaries. The fluid kinetic 
energy thereby depends on z. In the limit in which 
z >> a, the mass matrix ||Mjfc(z)|| of the sphere is well 
known [3; It is 

M + n ± {z) \ 

M + (i±(z) , (13) 

M + n { \(z) J 
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wherein 



M||W = m(i + |0 8 4 



(14) 



The potential energy of the sphere induced by fluid mo- 
mentum fluctuations is thereby 

U f i uid (z) = ^ lndet (m- 1 (z)M^) , (15) 
which reads 

3^k B T (a^ 3 

Ufluid(z) — — 
with // 



8(M + Vz 
27rpa 3 



and a « z. (16) 



From fluid velocity thermodynamic fluctuations, it fol- 
lows that the sphere will be attracted to wall with a 
potential proportional to the temperature and inversely 
proportional to the third power of the distance from the 
wall. 



wherein the statistical thermodynamic Hamaker 
constant CV is determined by the polarizability a? 



ksTar — d 2 — d 2 — d 2 (a « z), 

Ct — 4«t = 4a 3 (conducting sphere). (21) 

Altogether, the final attractive potential energy is 
fcgT /a\ 3 



U dipole (z) 



(-) +... (a«z). (22) 



Note the similarity between the fluid fluctuation po- 
tential in Eci. Ijl 6(1 and the dipole fluctuation potential 
in Ea.l(22|). Both potentials obey U oc -[k B T(a/z) 3 ] 
with proportionality constants of similar order unity. It 
would at this stage appear that the strength of electrical 
and fluid mechanical fluctuation forces are comparable in 
magnitude. However, this has only been proven at the 
classical statistical thermodynamic level of computation. 
Let us now consider quantum mechanical fluctuations. 



V. QUANTUM FLUCTUATIONS 



IV. ELECTRIC DIPOLE MOMENT 
FLUCTUATIONS 

If one places a neutral conducting sphere in the neigh- 
borhood of a perfectly conducting wall, then charge re- 
arrangements within the sphere will create fluctuating 
electric dipole moments. The fluctuating dipole moments 
will induce an attraction between the sphere and the wall 
as will now be shown. We again assume that the sphere 
radius is much less than the distance between the sphere 
and the wall a « z. The dipole moment d of the con- 
ducting sphere will induce a surface charge on the per- 
fectly conducting wall usually described in terms of an 
"image" dipole moment . The interaction between the 
dipole moment and the image is given by 



AF = 



d d, -3(d n)(d, n) 

ilk 3 ' 



(17) 



The image dipole moment is related to the dipole moment 
of the sphere via 

d z = d lz , d x = -d lx and d y = ~d iy , (18) 

so that the interaction free energy reads 



AF = 



4 



2d 2 



16z 3 



(19) 



Taking the thermodynamic average Udipoie — AF yields 



Udipoie \%) 



,+... (a«z) 



C T = d 2 x +d 2 y + 2d 2 z 



(20) 



In the quantum mechanical theory of fluctuations, the 
static response function k B TGij(Q) at zero temperature 
in Eq.lO is replaced by a complex frequency £ dependent 
response function k B TGij(Q, C) which obeys a dispersion 
relation with 3mC > of the form 

c ,, ft0 jr^i«*. (23) 



(« 2 - C 2 ) 

The static response function is then the zero frequency 
limit 

G tJ (Q) = hmG y -(Q,C) 

2 f 00 duj 
Gij(Q) = - SmG ij (Q,uj + iO + )— . (24) 

7T Jo u 

The power spectrum of quantum noise corresponding 
to the frequency dependent response function Gjj(Q,£) 
obeys the quantum mechanical fluctuation dissipation 
theorem 



1 



({AA^AA^O)}} 



S i3 (Q,u) e - luJt du, 



Sij(Q,u)) = ^^-^mG ij (Q,u; + iO + ), 



1 2 J \2k B T 



Employing an identity, 

E t {uj) 
k B T = 



(25) 



E 



or 



LU 2 + UJ 2 
n=—oo 

2ixk B Tn 



(26) 
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along with Eqs.J53, JUJ) and imply 



- (AXiAXj + AX, AX,) 



K B 1 



(27) 



n— — oo 



It is worthwhile to compare the quantum mechanical fluc- 
tuation Ea. (tT7|) to the classical Eq.Q; 



AXiAXj = k B TG t j(Q, 0) = k B TG %3 {Q). (28) 

For a single fluctuating variable, say X = a^A^, one 
obtains from Eq.J23 the expression 



(A* 2 ) = k B T E G(Q,ikn|), 



G(Q,C) = - 



n— — oo 
00 ,.,c>. 



jj MnGij(Q,uj + i0 + )du 



Employing an inequality, 
1 



< 



LO 2 -C 



1 



(29) 



(30) 



n n 



and a definition for the X— motion frequency scale ui a 
2 



into Ea. (|29|l yields the upper bound 

\(AX 2 ) - k B TG(Q,0)\ 

<2u; 2 oc k B TG(Q,0)Y / — , 



uj%mG(Q,uj + iO + )duj, (31) 



n— 1 '* 



i.e. 



(AX 2 ) 



k B TG(Q,0) 



< 



1 / 7iu; 



12 \ k B T 



(32) 



(33) 



From the inequality in Ea. lf^ we find that a sufficient 
condition for employing classical fluctuations is 



(34) 



(AX 2 ) « k B TG(Q, 0) if ^ « /c B T 
in agreement the classical Eq.©. 

VI. FLUCTUATION FREQUENCIES 



Our purpose is to estimate for both the fluid me- 
chanical and the electrical fluctuation forces. We con- 
clude for "room temperature" that 

TilUoq « k B T fluid mechanics (classical), 

hujoo >> k B T dipole moments (quantum). (35) 

The derivations follow. 



A. Effective Mass Sum Rule 

If p denotes the momentum of a colloidal particle 
within a fluid, then the dynamical mass of the particle 
obeys the Kubo formula 

M y (C) = -/ (IpMp^e^dt. (36) 



from which follows the equal time commutation sum rule 

— (\pi,Pj]) = - cj3mM„-(w + i0 + )duj. (38) 
n 1" Jo 

If the microscopic force p = f on the colloidal particle 
is derivable from a potential f = — gradV, then the sum 
rule EqEHI obeys 

2 f°° 

- / w3mM(w + i0 + )duj = (grad gradF) . (39) 
Ti" Jo 

On the other hand, from the dispersion relation, 



From Ea. l|36f) one finds 

l -{[ Pl {t), Pl m) 

> 

ZsmM i:j (uj + iO + )s'm(u;t)duj, (37) 



M(C) 



Lu^smM(u + iG + )dw 



7T JO ' C 2 ) 

it follows that the static mass obeys 



if 3mC > 0, (40) 



QmM(u + i0 + ) 



du) 



(41) 



From Eqs.JSni an d jlQ one computes the Hooks law fre- 
quency tensor 

fi L = M_1 ' (S rad gradV) . (42) 

For a given principle "X-direction" of the tensor, the hy- 
drodynamic frequency scale is given by 

2 1 /d 2 V\ 
00 M \ dX 2 J ' 



(43) 



The mass of the colloidal particle is proportional to 
the volume of the particle. The interaction V between 
the colloidal particle and the fluid is spatially non-zero 
only in the neighborhood of the particle surface. Thus, 
(d 2 V/ 'dX 2 ) is proportional to the contact surface area. 
If the number of atoms within the colloidal particle is 
denoted by N, then the number of atoms on the colloidal 
particles surface is proportional to N 2 / 3 . The frequency 
in Ea. l|43|l may be estimated by 



•jj- 



(44) 



where lo v h, is a typical atomic vibrational (say phonon) 
frequency. As a numerical example, let us consider a 
colloidal particle with N ~ 10 10 and with a vibrational 
frequency obeying hio v n, ~ k B T at room temperature. 
For such a colloidal particle hcu^ ~ 0.02 k B T which is in 
the classical regime of Eq. (13511 . 
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B. Polarizability Sum Rule 

If d denotes the electric dipole moment of a colloidal 
particle within a fluid, then the polarizability of the par- 
ticle obeys 

i r°° 

MC) = W WUjiO))) e«*dt. (45) 



From Eg . (|45|) one finds 



^smaij(uj + i0 + ) sin(u)t)du), (46) 

from which follows the equal time commutation sum rule 
2 



di, d 3 



u)$smotij(uj + iO + )du>. (47) 



The electric dipole moment and its rate of change, 
summed over all the charges within the colloidal particle, 
is given by 



ey^Zfci-fc, 

k 



(48) 



From the equal time commutation relation 
[vfc,rj] = —ihSki(l/Mk) one finds that [di,dj] = 



-ihSij Yjk^ z V M *) yielding 



uj\sma, 



j(u + iO + )du> — Sij 



From the dispersion relation (with 3to£ > 0) 
2 f°° tu^smaij(uj + z0 + )dw 



a«(C) 



(J* - C 2 ) 



it follows that the static polarizability obeys 

2 f°° du> 
otij = ay(0) = — ^smaij(u + i0 + ) — . 

7T Jo ^ 

For a spherical colloidal particle, 
2 



(49) 



(50) 



(51) 



Tra T J Jo 



uj^mu^ + «0 )dw. (52) 



For a conducting sphere of volume V — (47ra 3 /3) the 
frequency scale 



--T 



Aire 2 



E 



(53) 



The plasma frequency for a portion of condensed matter 



a 



2 _ ^ 2 S^fA. 



E 



ft 



(54) 



is dominated by electronic oscillations f2 2 « (47rne 2 /m) 
wherein to and n represent, respectively, the electron 
mass and density of electrons per unit volume. The fre- 
quency scale for dipolar fluctuations is then 



3 



(55) 



A metallic plasma frequency is of order (h^lp/kg) ~ 
10 5 °K. Eq.JSSJ) tnen implies hujoc >> k B T, as in 
Eq. (|35() . for temperatures near room temperature. 



VII. QUANTUM DIPOLAR FORCES 

It has been found at room temperature that fluid fluc- 
tuation forces are classical and electric fluctuation forces 
are quantum mechanical. An estimate of the quantum 
mechanical dipolar potential follows. The energy of in- 
teraction between a conducting sphere and wall due to 
dipole quantum fluctuations is found by summing the 
polarizability over Matsubara frequencies 8] as in SeclVl 



(56) 



A simple model for the polarizability will be found in 
order to estimate the potential energy in Eq.JSSJ. We as- 
sume that the polarizability has a single pole at frequency 
ujoc] i.e. 



"(0 



OiTUJr. 



c 2 



(57) 



The residue at the pole has been fixed so that a(0) = 
olt- Substituting Ea. (|57|l in Ea. H56fl one finds for the 
interaction potential 



Udipole{z) 



Udipole (^) 



k B T 

4z 3 

/ hujooar 



E 



V 8z 3 

Since fiw^ 3> k B T one finds 

Udipole\Z) — 



coth( 



2k B T 



,a 



)• (58) 



(59) 



By comparing Ea. (|59|l to Ea. 1)16(1 . it seems that the 
dipole fluctuation forces are much larger than the fluid 
fluctuation forces we find that 



Ufluid{z) 
U dipole (^) 



3/i 



M + fj, 



k B T 



(60) 



The term on the right hand side of Ea. (|60fl which is in 
square brackets is of order unity. Thus 



Ufluidjz) 
Udipole{ z ) 



kuT 



« 1. 



(61) 
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Using typical plasma frequencies for metals, Woo = 
(Sl p /y/3) ~ (10 16 /sec) or equivalently (/w^oo/fee) ~ 
10 5 °K. The inequality in Ea. l|61JI holds by a very large 
margin. 



VIII. FLUID MOTION 

The static classical fluid force between two spherical 
colloidal particles separated by r may be shown to be 
derived from the potential 



Ufi U id(r) = -37T 2 



k B Tp 2 a 6 
(A» + Mf 



(^) 6 (r»a), (62) 



wherein p is the fluid mass density and p, is the Eulcr 
mass and M is the bare mass of the colloidal particle. 
The above potential is derived from a long time scale 
statistical averaging over those fluid mechanical fluctua- 
tions which otherwise induce colloidal particle Brownian 
motion. Direct observation [10J of Brownian motion 
forces require shorter time scales. Typical experimen- 
tal bandwidths for micron scale colloidal particle size are 
about a tenth of a megahertz. 

For two identical metallic colloidal particles separated 
by r, the quantum electric field fluctuation potential is 
given by 



3/c T c *" ! 

Udipoie(r) = — ^ tt(i\oj n \)a(i\u> n \) 



71— — OO 



for r » a. 



(63) 



In the single pole approximation for a(£) with the metal- 
lic a(0) = a 3 we find that 



V r 



Udipole(r) ~ - 

for r » a and Tilj^ >> k B T. 

Note that 

Ufiuidjr) 



Udipole(r) 



p 2 a 6 



{n + My 



k B T 



flUJr. 



« 1. 



(64) 



(65) 



The very large magnitude of ~ (10 16 /sec) forbids 
direct observation of dipolc fluctuations. 



IX. CONCLUSIONS 

The general theory of quantum mechanical fluctuations 
was discussed with particular emphasis on computing the 
frequency scales of motion from sum rules. The frequency 
scales determine whether the fluctuations and thereby 
the forces are classical or quantum mechanical in nature. 

We have shown how fluctuations in fluid velocity fields 
and in electric fields give rise to forces exerted on colloidal 
particles. Fluctuation forces were computed in detail for 
the case of colloidal particles attracted to the walls of the 
suspension container. The resulting van der Waals force 
has the form 



Uwaii(z) = -Ui (a/zf for z >> a. 



(66) 



The electric field fluctuation contribution to U± dom- 
inates the fluid mechanical contribution to U± as in 
Ea. (|61|) . The long ranged static attraction between two 
spherical particles of radius a separated by a distance r 
has the van der Waals form 



U,, 



rticl 



,(r) = —Uq (a/r) for r >> a. 



(67) 



The electric field fluctuation contribution to Uo dom- 
inates the fluid mechanical contribution to Uq as in 
Eq.©. 

At the level of classical statistical thermodynamics, the 
fluid velocity and electric field contributions to this po- 
tential are comparable. When quantum fluctuation ef- 
fects are taken into account, the electric fluctuation con- 
tribution to the potentials dominates the fluid mechani- 
cal contribution to the potentials as in Eas. l|61|l and 16511 . 
The electric field fluctuation long ranged forces can be, 
and have been|n],0], observed by measuring phase sep- 
arations in some colloidal suspensions. Further work on 
colloidal forces in still other geometries would be of gen- 
eral interest. 
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